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Abstract 

We illustrate the equivalence between the non-unitary evolution of an open quantum system governed 
by a Markovian master equation and a process of continuous measurements involving this system. We 
investigate a system of two coupled modes, only one of them interacting with external degrees of freedom, 
represented, in the first case, by a finite number of harmonic oscillators, and, in the second, by a sequence 
of atoms where each one interacts with a single mode during a limited time. Two distinct regimes appear, 
one of them corresponding to a Zeno-like behavior in the limit of large dissipation. 
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The decoherence program is one of the most successfull theroretical proposals devised to give 
suitable responses to the problem of quantum-to-classical transition (see |ID and references therein) 
within the scope of the quantum theory itself. According to this program, isolated systems must 
be considered as mere idealization and the emergence of classicality from quantum mechanics 
is a consequence of the unavoidable interaction between a physical system and its environment. 
The net effect of this interaction is the disappearance of the local quantum correlations between 
preferred states^. The early formulations of this program assign the mechanism of decoherence 
to the continuous monitoring of the system of interest by the surrounding particles |2j-l5|]. These 
formulations were devised to explain the quasilocalization of macroscopic bodies or symmetry 
breaking in microscopic systems, like the observation of a well defined chirality in optical isomers. 

One of the widespread ways to obtain the effective dynamics of an open quantum system con- 
siders the microscopic model describing the coupling between the system of interest and the de- 
grees of freedom of its environment. In general, the number of these environmental degrees of 
freedom considered is very large, such that the irreversibility of the effective dynamics could be 
guaranteed. In the limit of infinite external degrees of freedom, the environment is suitably mod- 
eled by a thermal reservoir and the effective dynamics of the system of interest is governed by a 
master equation obtained by tracing these external degrees of freedom out flg]. 

Another well known effect related to continuous measurements is the Quantum Zeno Effect 
(QZE). In 1977 B. Misra and E. C. G. Sudarshan reported an intriguing result related to mea- 
surements in Quantum Mechanics |7|]. They showed that a sequence of projective measurements 
on a system inhibits its time evolution. In the limit of continuous measurements the evolution is 
completely frozen. Similarities with one of the paradoxes proposed by Zeno of Elea, who intend 
to show that movement is theoretically impossible, motivated the authors to name the quantum 
effect after the Greek philosopher. Originally the Quantum Zeno Effect was called the Quantum 
Zeno Paradox, because it was supposed to show the theoretical impossibility of a "movement" (the 
quantum evolution) of a decaying particle in a bubble chamber. 

The Quantum Zeno Effect (QZE) plays an important role in Quantum Mechanics and enormous 
literature about this topic has been produced over the last thirty years. After the realization of the 
pioneer experiment [|8D on the effect, which showed the inhibition of transitions between quantum 
states by means of frequent observations, the QZE became the center of fervorous debates ^.[hI 
1 To be exact, the quantum correlations does not disappear. In fact, they are carried from the system to the environ- 
ment and turned out inaccessible. 
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about the necessity of t 
have been proposed [9, 



le projection postulate on the measurement description. New approaches 



1 IN and the initial association between the QZE and the projection postu- 



late was no longer a necessary ingredient. Nowadays the literature on this subject is extensive and 



includes relation between QZE and quantum jumps |. 
implementation in the system of microwave cavities 



130, quantum Zeno dynamics II 14 , 



150. its 



161, and semiclassical evolution for coupled 



systems obtained by frequent Zeno-like measurements II 1711 . With the increasing interest in quan- 
tum computation, QZE has become also a tool for the development of protocols on quantum state 



protection Ill8l-l20ll. that are important for the implementation of quantum computation. 

Although originally the QZE was considered as a result of a sequence of measurements on the 
system of interest, the effect has also been predicted in different contexts, as an example we quote: 
In Ref. n2lh the authors suggested that a manifestation of the QZE appears in the relaxation of 
optical activity in a medium. In Ref.|220 a relation between the QZE and the coherence properties 



24|] the 



of a two-state system subjected to random influences in a medium is presented. In Ref.1123 
authors also show that the QZE can be induced by other physical interactions. 

In the present work we study two completely distinct approaches to the dynamics of quantum 
open systems. We consider the system composed by two coupled modes interacting unitarily, 
where only one of them is coupled to external degrees of freedom. We show that the Zeno like 
effect similar to the one shown in ||24l is present in both dynamics. Firstly we couple to the system 
of interest a dynamical environment over which we have control (interactions, number of com- 
ponents, etc.). Secondly we simulate an experimental situation where the system is continuously 
monitored by a probe system, again controlling interaction time, and other parameters. 

We show that not only the conjecture is correct in the situations presented, but also a Zeno 
like effect can be obtained for both cases, when the influence of the external system is sufficiently 
strong. Precisely the same effect has been observed by modeling the tunneling of a photon between 
two cavities, one of them whose dissipation is governed by a master equation [|2i 



L COUPLING WITH A FINITE NUMBER OF HARMONIC OSCILLATORS 

In this section we consider the system of two linearly coupled harmonic oscillators, one of them 
coupled to an environment composed by a finite number of harmonic oscillators according to the 
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Hamiltonian 



H = H s + H R + Hi nt , H s = coia[di + u\a} 2 a 2 + G ( a\a 2 + a\hi 



(1) 



N N 

H R = ^2oj{a[a k , H int = (a\ 



fc=3 fc=3 

where a\ and a\ (Si and a?) are creation (annihilation) bosonic operators for the modes of interest 
Mi and M 2 , and a k and a k refer to the environment for k ranging from k = 3 to N. Defining 
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we can write the Hamiltonian in a matrix form: 



# = A^HA. 



(3) 



In order to preserve the hermiticity of the operator H, the matrix H must be hermitian. Without 
loss of generality, we can consider it real and symmetric. Thus, there is an orthogonal matrix P 
such that 



P 7 HP 



^ Ai ••■ ^ 
A 2 •■■ 



y • • • X N j 



D, 



(4) 



where the \ k are real numbers, which may be used to write the Hamiltonian in a diagonal form: 



M 



(5) 



k=l 



where 



B = P A 



B D = P T A D 



(6) 
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Using the orthogonality of P, it is easy to show that canonical commutation relations hold for b\ 
and b\: 



h,b\ 



Skj, 



b k ,bj 



0. 



(7) 



The modes related to the operators shall be called the original modes; the ones concerning to 
b k will be referred to as the normal modes. 

In what follows, we assume that the main modes Mi and M 2 are resonant and we consider that 
the matrix H is a function of five positive parameters N, f2, G, A and T: 
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The parameter f2 gives the frequency of the oscillators of interest, which are coupled accord- 
ing to the constant G. The total number of environmental oscillators is (N — 2). Their fre- 
quencies are distributed in a finite interval around Q defined with the help of the parameter A: 
[Cl — A/2, f2 + A/2]. If A is large enough, the environmental modes out of this interval have 
negligible effects over the oscillators of interest [250. The environmental frequencies are equally 
spaced distributed above and below the central frequency f2. This choice is not essential for the 
dynamics induced and for derivation of master equations: in fact, randomly distributed frequen- 
cies lead, when N is large, to equivalent results. In derivations of master equations, it is usually 
assumed reservoirs with infinite frequency range and dense spectrum. Infinite frequency range is 
approximated by increasing A, and dense spectrum by decreasing A/ (N — 2); in order to achieve 
both conditions simultaneously, we must have large values of A and N. In this limit, a master 
equation may derived from the Hamiltonian specified by Eq. ©. The parameter T is related to the 
strength of the coupling with the environment. The choice of 7^ being proportional to 1/ y/N — 2 
is consistent with assumptions usually performed in derivation of master equations that guarantee 
finite decay rates in the thermodynamic limit N — > 00. Of course, all 7^ do not have to be equal 



for such derivations. 

Since the b\ are linear combinations of the fife, they share the same vacuum state |0). In order 
to investigate the dynamics of the system plus environment in the space of one excitation, let us 
define the states 

|0 fc ) = 4l°>. \<Pk) = b{\0), k = ltoN. (9) 
Using Eqs. ©, we see that they are connected by 



N N 



(10) 



where Pjj is the element of the matrix P in the z-th line and j-th column. Using these relations and 
observing that H \ <j> k ) = X k \4>k), we can calculate the evolution of the states with one excitation 
in the original modes as 

N N / N \ 

e~ lAt \e k ) = J2 p ^ e ~ im lfc> = E H p ^ e ~ lXjt W ■ ( n ) 

j=l 1=1 \j=l J 

The probability of finding one excitation in mode M 2 if it is initially in this mode is 



P(t) 



N 



A, I 



(12) 



In Ref. 112411 . two regimes were found for the dynamics of such a probability calculated by using 



the master equation corresponding to Hamiltonian (OQ) with resonant modes M 1 and M 2 : 

—p = k (^2aipa\ — pa\ai — ajdipj (13) 



tin 



a\di + dg^, p 
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d 2 ai + a\a.2, p 



Here, the density operator p stands for the state of the composed system Mi + M 2 . In one of 
the regimes, the increasing of the dissipation constant n of mode Mi leads to the decreasing of 
the permanence of the excitation in mode M 2 . This is expected, since mode M 2 is connected to 
the environment only through mode M±. The other regime was called Zeno regime: there, the 
increasing of the interaction of Mi with the environment inhibits the transition of the excitation 
from M 2 to Mi, leading to the enhancement of the probability of finding the excitation in M 2 . 
The turning point between these regimes occurs for k = 2G, where G is the unitary coupling 
constant between the modes of interest. In the Appendix, we show that such a turning point 
corresponds to T = y / 2AG/V, what is corroborated by Figs. (fTal) and dial) . The occurrence of 



two regimes may be understood with the help of the following analytically calculable cases: for 
T = 0, p (t) = (1 + cos (2Gt)) /2; for G = 0, p (t) = 1. Depending on the relation between 
T and G, the dynamics can be approached to one of these limiting cases. By comparing Figs. 
(fTal ) and (|2ab with Figs, (flbl) and (|2bl , respectively, we see that the Hamiltonian and the master 
equation approaches exhibit good agreement. In order to achieve such an agreement, we have 
to pay attention to two aspects concerning the parameter A: it must be large enough so that 
environmental modes with frequencies out of the interval [Q — A/2,fi + A/2] have negligible 
action on the system (as pointed out in Ref. [25], the relevant environmental modes are the ones 
with frequencies around the frequencies of the normal modes of the system f2 ± G); the ratio 
A/ (N — 2) must be small, allowing the use of the limit of dense spectrum. 



II. SEQUENCE OF MEASUREMENTS 

In this section we study the dynamics of two resonant coupled modes (Mi, M 2 ) and N atoms 
interacting, one at the time, with mode M\. The sequence of interacting atoms represents, in the 
limit of instantaneous interactions, a continuous measurement of the M 1 excitation number. The 
investigation shows that the two regimes reported in Ref. J24I and in the previous section are also 
present if a sequence of atomic interactions is performed on mode M±. These results illustrate the 
relation between continuous measurement of a quantum system and the dynamics governed by a 
master equation. 



The results are obtained by numerical simulations, where we consider, as in Ref. [24], the 
system of modes Mi and M 2 in the initial state |0i, 12). A sequence of two level atoms interacts 
with mode M\. The atoms are prepared in the ground state \g) and interact with mode Mi one at 
the time. The Hamiltonian of the global system for the interaction of the A;-th atom is given by 

= n(a\a x + a\d 2 ) + G{a\a 2 + ajai) + ~<jW + g{d^0 + ai<rf } ), (14) 

where a\ (ai) and a\ (S2) are creation (annihilation) operators for modes M\ and M 2 , Q their 
frequency, G the modes coupling constant, 5T ' = |e^)(e^| — Ig^Xg^l, cri*^ = Ig^Ke^l, 
<3"+ = \e^) (g^ I and g the coupling constant for the interaction between k-th atom and mode M\. 
Here, \g^) and |e^ fc ^) stand for the ground and the excited states of the A>th atom, respectively. 

After each interaction we perform the trace over the corresponding atomic system, i.e., we do 
not consider the final state of the atoms. We also assume that the coupling constant g scales as 
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1 / y/ti n t, where t int is the interaction time of each atom [26]. The overall effect of these atomic 
interactions with a cavity mode is a dissipative evolution of the mode. The effective dissipative 
constant related to this process is k = g 2 t int . 

In Figs. (ITcb and (|2cl we shown two regimes for the probability of finding the excitation in 
mode M 2 (V(t)). In the dissipative regime (k < 2G), the increasing of the dissipative constant k 
leads to the decreasing of V(t). In the second regime (k > 2G), the increasing of the dissipative 
constant k leads to the increasing of V(t), preserving then the excitation in mode M 2 . It is worth 
to note that the agreement with the master equation results is reached in the limit of vanishing 
interaction time, t int . The increasing of interaction time leads the curves away from the ones 
obtained by the master equation. 



HI. CONCLUSION 

In the present work we investigate the two regimes of the system of two coupled modes, in- 
duced by two completely different dynamics. In the first one we consider the mode Mi linearly 
coupled to a finite number of harmonic oscillators, and in the second one we consider such mode 
interacting with N atoms, one at the time. Both dynamics, in appropriate limits, can describe the 



regimes obtained in [24] using a master equation. In the first model, when the number of harmonic 
oscillators goes to infinity the coupling between them and the system of interest can simulate the 
interaction with the environment. In the second model, the interaction with N atoms, when the 
interaction time goes to zero, simulates continuous measurements. Therefore, the present results 
illustrate the idea that the interaction between system of interest and environment can be inter- 
preted as continuous measurement on the system of interest. The results for both dynamics were 
obtained using numerical simulations. As no approximations were used in the calculations, the 
analysis of intermediate scenarios, out of the master equations limits, is possible. 
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Appendix A: Establishing the relation between the effective dissipation constant and the parame- 
ters used in the Hamiltonian approach 



The master equation employed in Ref. 1124(1 is given by Eq. (fT3l) . It describes the dynamics 
of two linearly coupled resonant modes, one of them interacting with an environment at zero 
absolute temperature. Such a master equation may be derived as an approximation of the dynamics 
emerging from the Hamiltonian H, leading, under the specifications in Eq. ([8]), to 
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where Uj = jA/ (N — 2) and r c is a value of |t| beyond which the summations above are negligi- 
ble Il25n . For N sufficiently large, we can take the limit of dense spectrum, resulting in 
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Since the integrand is assumed to be negligible for |r| > r c , we change ±r c for ±oo, which leads 
to 
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Observing that the term related to 5 (u + G) vanishes, we find, for A/ (N — 2) < G < A/2, 
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The transition to the Zeno regime were found in Ref. [24] for k = 2G. This corresponds to 
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(a) 



(b) 



(c) 



FIG. 1 : (a) Evolution of the occupation of mode M 2 according to the Hamiltonian approach for 
N = 500, Q = 100G, A = 20G and T = J^f- with 77 = 2 (solid line), 77 = 8 (dashed line) and 

77 = 32 (dotted line), (b) Evolution of the occupation of mode M 2 according to the master 
equation approach for k/G = 2 (solid line), k/G = 8 (dashed line) and k/G = 32 (dotted line). 

Figure from Ref. j24 \. (c) Evolution of the occupation of mode M 2 during a sequence of 



measurements probing the presence of the excitation in mode Mi for g — \ f 2 ^ with 77 = 2 



(solid line), r\ = 8 (dashed line) and 77 = 32 (dotted line). For sake of comparison, in the plots (a) 
and (c), a curve concerning a given value of rj corresponds to the curve related to the same value 

of k/G in plot (b). 
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(a) 



(b) 




(O 



FIG. 2: (a) Evolution of the occupation of mode M 2 according to the Hamiltonian approach for 



N = 500, n = 100G, A = 20G and T = J^- with 77 = 2 (solid line), 77 = 1/2 (dashed line) 
and 7] = 1/8 (dotted line), (b) Evolution of the occupation of mode M 2 according to the master 
equation approach for k/G = 2 (solid line), k/G = 1/2 (dashed line) and k/G = 1/8 (dotted 
line). Figure from Ref. Il24ll . (c) Evolution of the occupation of mode M 2 during a sequence of 
measurements probing the presence of the excitation in mode Mi for g = \Jjf^ with 77 = 2 
(solid line), 77 = 1/2 (dashed line) and 77 = 1/8 (dotted line). For sake of comparison, in the plots 
(a) and (c), a curve concerning a given value of r] corresponds to the curve related to the same 

value of k/G in plot (b). 
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